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Abstract 

The evolution of a two-state system driven by a sequence of imperfect vr 
pulses (with random phase or amplitude errors) is calculated. The resulting 
decreased fidelity is used to derive a plausible limit on the performance of 
"bang-bang" control methods for the suppression of decoherence. 
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I. INTRODUCTION AND MODEL 



A number of dynamical methods for the suppression of decoherence in quantum systems 
have recently been proposed . While some of these are quite sophisticated, the simplest 
such approach, sometimes called "bang-bang" control is a rather straightforward idea 
which could, in principle, be applied to any two-state system interacting with an environment 
which has a finite (nonzero) correlation time r^. (See ^ for a recent proof-of-principle 
experimental demonstration.) The idea is to "flip" rapidly (faster than Tc) the state of the 
system, back and forth, in such a way that the environment's unwanted influence on the 
system is constantly being undone by the environment itself. 

For instance, for a system to environment coupling which involves only the operator 
(pure phase decoherence), due to a term such as 

fi(^z (gkhi + alh) (1) 

k 

in the Hamiltonian (where the boson operators hk represent modes of the environment), a 
sequence of very short vr pulses, which rotate the system's pseudospin by 180° around the x 
axis and hence change the sign of cr^, would result in an evolution in which the sign of (|l]) 
changes form one instant to the next, and its effect therefore averages to zero. 

The purpose of this note is to consider the constraints which imperfections in the vr 
pulses place on the successful implementation of such a dynamical decoupling strategy. The 
approach is simply to consider the evolution of a two-state system subject to a series of 
imperfect pi pulses, and to calculate the total error (as measured by the state's fidelity) 
introduced, on the average, after such pulses. The average considered here is over all 
possible initial states of the system, as well as over the distribution of the random errors in 
the TT pulses. 

I work in the interaction picture and assume that the central frequency of the pulses 
is exactly tuned to the resonance frequency of the transition between the two states of 
the system. The rotating-wave approximation is also assumed to apply. The Hamiltonian 
describing the interaction of the system with a single pulse is then 

H = ng{£{t)a^ + £*{t)a^) (2) 

where 0"+ = |+)(— |, cr_ = |— )(+|, and |-|-) and |— ) are the eigenstates of a^. The coupling 
constant g can be taken to be real without loss of generality. In general, the Hamiltonian 
(0) does not commute with itself at different times. To simplify matters, I shall take the 
classical field amplitude to be of the form £{t) = E{t)e^'^, where E{t) is real and is a 
constant, determined essentially by the timing of the pulse. Then, the evolution operator 
corresponding to (H) is, in the {|+),|— )} basis. 



U{t)= exp 
= exp 



H{t')dt' 



h Jo 

—i (cos (fax — sin Lpay) J gE{t')dt' (3) 



If the duration of the pulse is tp, define the angle 6 = 2 /q'' gE{t')dt'. In pseudospin terms, 
the operator U (tp) represents a rotation of the spin by an angle about an axis in the x — y 
plane which makes an angle ip with the x axis. If O = vr, one has a vr pulse. 
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In this paper I consider the state resulting from a sequence of "imperfect" vr pulses. 
Imperfections may arise from either the phase (p, which controls the axis of the rotation, and 
which would be sensitive to errors in the timing of the pulses, or from errors in the amplitude 
and/or duration of the pulse, which would make 9 ^ n. Amplitude errors are considered 
in Section 2, phase errors in Section 3, and the results, and their possible consequences for 
"bang-bang" control, are discussed in Section 4. 



II. AMPLITUDE ERRORS 

In this Section I shall assume that = 0, in Eq. (|^), for all the pulses in the sequence. 
A full spin cycle consists of two vr pulses (separated by a short time of "free" evolution), 
so I shall consider a total of 2A^ pulses, where is the total number of cycles. I will take 
the "pulse area" to be of the form O = vr + Cj, where is a small number which varies 
randomly from one pulse to the next. All the are assumed to be distributed in a Gaussian 
way, that is, according to the following probability distribution: 



/27rA 



(4) 



with standard deviation A. 

With this notation the evolution operator for the i-th pulse is simply 



(5) 



and the result of 2A^ consecutive pulses is given by 



total ~ 



-i{Ar7r+ei/2+...e2Jv/2)o-^ 



f-1 



e . . e 

cos I sm -cr^ 

2 2 



(6) 



where the random variable e = ei + . . . t2N is also distributed in a Gaussian way, and has 
standard deviation \/2NA. 

A general initial state of the two-level system can be written as 



9 9 

IV-o) = cos-|+) + e*^sin-|-) 



(7) 



In this state the spin is aligned along an axis of colatitude 9 and azimuth (p. The effect of 
(D on the state (0) is 



Utotalli^o) = (-1) 



N 
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9' 



cos - COS ie*"*' sin - sin - I +) + cos - sin -e^'^ — i sin - cos - I — ) 

2 2 2 2;''V2 2 2 2/'^ 

(8) 



The fidelity, JF, of the state after all the pulses is given by JF = \ {il^o\Utotai\'4^Q)^ , and it equals 



J- = COS — h sm - sm 9 cos i 



(9) 
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To calculate an average fidelity, Eq. (|]) may be integrated over 6 and 0, assuming a uni- 
form distribution of initial states over the whole Bloch sphere (that is, with a density 
sin 9d9d(l)/ in), and over e with the appropriate probability distribution. Before proceed- 
ing to do so, however, note that sin 6' cos is just x in Cartesian coordinates; hence a 
reparametrization of the Bloch sphere that exchanges x and z will turn sin 6* cos into cos 6', 
which makes the integral over (p' trivial. The probability distribution for the final state 
fidelity can thus be written formally as 

P(^)= -\^L= r e-^'/^^^'tie r sine'dd'5 ( T - cos^ - - cos^ Q' sin^ - 



VWn J- 

11 °° rVT 




"^"^ (10) 



{l-T){T-x^) 



The last integral may be evaluated numerically for a given value of A^A^ and T (with 
< JF < 1); for reasonable values of iVA^, only a few terms in the sum over ra, around n = 0, 
are necessary. The results are shown in Figure 1 for A^A^ = 0.1, 1, and 10, respectively. 
The divergence as — > 1 is integrable; the value of P{T) as — > is finite. The average 
of JF can be calculated easily, and the result is 

{:f) = 1 + ^e-^^^ (11) 

This shows the expected decrease in fidelity with the number of cycles, but it also shows that 
for sufficiently large iVA the average fidelity saturates at the value 2/3. This may look like 
a relatively large number, but one must keep in mind that for a two-level system a totally 
random mixed state with density operator p = 1/2 would still yield a fidelity of 0.5 when 
compared to any initial state. Thus, in fact, a fidelity of 0.67 for this system corresponds to 
an almost random final state. 

This is illustrated in Figure 2, which shows three simulations of trajectories starting 
from three randomly picked states and going through a total of N^ax = 400 cycles, with 
A chosen in each case so that N^ax^^ = 0.1, 1, and 10, as in Fig. 1 (the corresponding 
average fidelities, as given by Eq. (pH]) , are 0.968, 0.789 and 0.667). Clearly, the trajectory 
for \/Nrnax^ = 10 is all over the place. (The trajectories shown have been chosen from 
a relatively large sample of randomly generated trajectories so as to be as "average" as 
possible, that is, neither unusually "good" nor unusually "bad" for the particular value of 
Njnax^^ that they illustrate.) 

In this connection it may be worthwhile to consider the "worst-case fidelity" predicted 
by Eq. (|^): it corresponds to initial states having either ^ = 0, vr, or = 7r/2, and, when 
averaged over e, yields 

{^rmn) = \ + ^e"^^^ (12) 

For large NA"^, this gives 1/2, i.e., the fidelity of the totally random mixed state, as discussed 
above. 
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III. PHASE ERRORS 



In this Section I take the pulse area G = 2 Jq^ gE(t')dt' to be exactly vr, and the phase 
if in Eq. @ to be a small random number, varying from pulse to pulse, with zero average 
and a Gaussian distribution with standard deviation A. Then the evolution operator for 
the i-th pulse is 

JJ.— ^-i{TT/2){cOSipiaa:SinipiUy) 

= — z(cos vjjCXa; — simpiay) (13) 
and its effect on the general state (|^ is 

f/#o) = e-*^' (cos ^1+) + e-'('^+2'^') ^1")) (1^) 

that is, except for an overall phase factor, all that happens is that the angle becomes 
(j) + (fi. Therefore, after cycles totalling 2N pulses, the final state will be, up to a phase 
factor 

Utotaim = cos^|+) + e*('^+^)sin^|-) (15) 

where e = 2{lpi + . . . + ip2N) is now a random variable with zero average and a Gaussian 
distribution with standard deviation 2\/2NA. A little trigonometric manipulation shows 
that the fidelity of the state ([T5|) can be written as 

jr = cos^ - + sin^ - cos^ 9 (16) 

which is of the same form as (^ (after the reparametrization of the Bloch sphere which 
exchanges x and z), and yields the same results (|r^) and ( |TI| ) for the probability distribution 
and average value of the fidelity, only with replaced by 4N everywhere. 



IV. DISCUSSION 

From the analysis in the previous sections one can conclude that, if each pulse used for 
"bang-bang" control has a r.m.s. error of A radians in either the phase or the (integrated) 
amplitude, in order to avoid loss of fidelity one must (roughly speaking) restrict the total 
number of cycles to 

Nmax < (17) 

In practice, of course, "bang-bang" control would be applied to a system which has an 
internal Hamiltonian and is, as well, interacting with other systems, so the overall picture 
will be rather more complicated than that described by Eq. (H), but there are grounds to 
believe that the conclusion above should still hold true, approximately, in spite of these 
complications. For instance, recent calculations of decoherence due to internal interactions 
in a quantum computer show that the net effect does not depend much on whether the 
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computer is actually carrying out a calculation or merely sitting idle in a predetermined 
state H]. Similarly, Miquel et al. [0 found, in their simulation of a quantum computation 
with phase drift errors, that the effect of said errors could be well approximated by a formula 
which did not depend on the details of the computation, only on the number of pulses. 

As a simple test, I show in Fig. 3 the results of a simulation for a system with the 
self-Hamiltonian 

Ho = hQa, (18) 

initially prepared in the state {ipo) = 75(1+) + which, without "bang-bang" con- 

trol, would evolve in time as {ipt) = + ie*™|— )), and hence exhibit a fidelity 

= \{ipo\Utotai\'ipo)f = cos^ Qt (solid curve). Low-noise "bang-bang" control (dashed curve, 
amplitude noise, Nmax^"^ = 0.1) cancels out this free evolution quite effectively, so the state 
fidelity remains close to 1 for the times shown, but if the noise is high instead (dotted curve, 
Nmax^"^ = 10), the fidelity is completely lost. This simulation assumes that the control 
pulses are essentially instantaneous and very strong, so that the self-Hamiltonian ([T8|) can 
be ignored during a control pulse. Also, in order to keep the fidelity close to 1, the interval 
between pulses At needs to be short enough; specifically, in this case, shorter than 2Tr/Q. 
For Figure 3, At = O.OOSvr/fi. 

It should be pointed out that the initial state chosen for the simulation in Fig. 3 (an 
eigenstate of ay) is just about the "worst" possible choice, that is, it is one of the states 
that minimize the fidelity as given by Eq. (^ (the average fidelity for such states is given 
by Eq. ([I^)). This explains why the simulations in Fig. 3 look somewhat worse than those 
in Fig. 2. Nonetheless, this initial state was chosen because it gives maximum modulation 
of the fidelity under the self-Hamiltonian ([T8|) . 

In Ref. it was shown that, in order to combat decoherence by the "bang-bang" method, 
one needs to have at least one complete cycle every Tc, where Tc is the characteristic "coher- 
ence time" of the environment. Putting this condition together with (|T^ yields the result 
that the "bang-bang" method with imperfect pulses can protect a quantum information 
storing (or processing) device for a time no longer than 

where A is the r.m.s. error per pulse. For instance, if the vr pulses are known to be accurate 
to, say, one part in 10^, one may expect A ~ tt x 10^^, and Tmax ~ IO^Tc 

It is important to keep in mind that Tc is not, in general, equal to the decoherence time 
of the system itself under its interaction with the environment. As an extreme example, 
the spontaneous emission of an atom has a decoherence time of the order of I/7, where 7 
may be mega- or gigahertz for an optical transition, but the environment responsible for the 
decoherence — namely, the quantized electromagnetic field — has a "coherence time" of the 
order of an optical period or shorter, that is, about 10~^^ s. On the other hand, for other 
kinds of decoherence — such as, for instance, the kinds resulting from internal interactions 
between the qubits — one may expect the relevant Tc to be of the order of a characteristic 
interaction frequency, which would also determine the basic decoherence timescale. In that 
case, "bang-bang" control with accurate pulses could substantially lengthen the decoherence 
time, as Eq. ([ToD indicates. 
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FIGURES 



FIG. 1. Probability distribution of the fidelity for a totally random (uniformly distributed on 
the Bloch sphere) initial state, and Gaussian noise of standard deviation A per pulse, for different 
values of iVA^. SoHd hue: TVA^ = 0.1. Dashed line: iVA^ = 1. Dotted line: iVA^ = 10 

FIG. 2. Fidelity as a function of number of pulses for different values of A. Solid line: 
Nmax^^ = 0.1. Dashed line: Nmax^^ = 1- Dotted line: Nmax^^ = 10. In ah cases Nmax = 400. 



FIG. 3. Fidelity as a function of number of pulses for a system with the self-Hamiltonian (18) 
for different values of A. Solid line: free evolution, no pulses. Dashed line: N^ax^^ = 0.1. Dotted 
line: N^ax^'^ = 10. In all cases Nmax = 400. 
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